1 



o 
o 
o 

(N 
> 

o 



> 

On 
(N 
O 



O 

o 

-f— > 
Oh! 



/b f° r Various Actions : Approaching the Continuum Limit with 
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We present results for pseudoscalar decay constants of heavy-light mesons using both quenched and Nf = 2 
dynamical fermion configurations. A variety of fermion actions is investigated: Wilson, nonperturbative clover, 
and fat-link clover. For heavy quarks the Fermilab formalism is applied. In the quenched approximation, results 
with the nonperturbatively improved clover action of the Alpha collaboration allow us to study the systematic 
error of the continuum extrapolation from the Wilson action. In addition, we use quenched configurations to 
explore the effects of fattening. The lessons from the quenched analyses are then applied to data with dynamical 
fermions, where both Wilson and fat-link clover actions have been used. This allows us to attempt a continuum 
extrapolation of the dynamical results. 
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Decay constants of the B and B s mesons are 
essential for accurate determination of parame- 
ters of the CKM matrix. We present here lat- 
tice studies of decay constants of heavy-light 
mesons. A MILC collaboration computation of 
the quenched decay constants has appeared pre- 
viously jjj; in that work, dynamical simulations 
were used only to estimate the quenching error. 
Our first continuum-limit results for dynamical 
lattices with Nf = 2 staggered quarks were pre- 
sented at Lattice '99 [|| . Here the dynamical lat- 
tices are studied in more detail, with both Wil- 
son and clover valence fermions. We also update 
the quenched lattice results, based on some new 
runs and revised analysis. Some of these new runs 
use clover fermions with different levels of fatten- 
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ing (including no fattening, the standard "thin" 
case), thereby enabling us to understand better 
both the fat-link clover dynamical results and the 
quenched discretization errors. 

Table @ gives the lattice parameters. The sets 
J and CP are new additions to the quenched lat- 
tices previously analyzed in Ref . . Unimproved 
Wilson valence fermion propagators were gener- 
ated for all these sets except J. On J, and a 200- 
lattice subset "CP1" of CP, we used thin link 
clover fermion propagators. For a 100-lattice sub- 
set "CP2" of CP1, we generated clover propaga- 
tors with four different levels of fattening (see be- 
low). For the dynamical lattices, valence Wilson 
fermions were studied on all the lattices; fat link 
clover fermions were also studied on a 98-lattice 
subset "RF" of set R, with 10 levels of fattening 
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and coefficient c = 0.45. 

For our central values, the dynamical quark 
configurations are treated as fixed backgrounds 
and chiral extrapolation is performed in the va- 
lence quark mass only. To estimate the system- 
atic error due to this "partial quenching," we 
also extrapolate with m va i e nco = ^dynamical, the 
equality being defined by parameters correspond- 
ing to the same m v . 

Table 1 

Lattice parameters. The upper group corre- 
sponds to quenched lattices; the lower group, to 
dynamical lattices with Nf = 2 staggered quarks. 
The set G was generated by HEMCGC. 



name 


P(am g ) 


size 


# configs. 


A 


5.7 


8 3 x 48 


200 


B 


5.7 


16 3 x 48 


100 


E 


5.85 


12 3 x 48 


100 


C 


6.0 


16 3 x 48 


100 


CP 


6.0 


16 3 x 48 


305 


J 


6.15 


16 3 x 48 


200 


D 


6.3 


24 3 x 80 


100 


H 


6.52 


32 3 x 100 


60 


L 


5.445 (0.025) 


16 3 x 48 


100 


N 


5.5 (0.1) 


24 3 x 64 


100 





5.5 (0.05) 


24 3 x 64 


100 


M 


5.5 (0.025) 


20 3 x 64 


199 


P 


5.5 (0.0125) 


20 3 x 64 


199 


G 


5.6 (0.01) 


16 3 x 32 


200 


R 


5.6 (0.01) 


24 3 x 64 


200 


S 


5.6 (0.02) 


24 3 x 64 


202 


T 


5.6 (0.04) 


24 3 x 64 


201 


U 


5.6 (0.08) 


24 3 x 64 


202 



For clover fermions, both thin and fat, the full 
Fermilab formalis m [p| is applied, including use of 
the EKM norm yl — 6k, the 3-d rotation of the 
fields, and the identification of kinetic mass mi 
(rather than the pole mass mi) as the physical 
mass. The shift m\ — > ni2 is performed at tad- 
pole improved tree level. In the thin-link case, we 
take the clover coefficient csw and the light-light 
renormalization/improvement constants from the 
nonperturbative determinations by the Alpha col- 
laboration For the heavy-light renormaliza- 
tions around the B meson mass, we use either the 
one-loop perturbative calculations H or a ver- 



sion of tree-level tadpole improvement designed 
to smoothly join on to the nonperturbative results 
as the quark mass gets small. The former ap- 
proach is designated "NP-IOY;" the latter, "NP- 
tad." For more details on NP-tad, see Ref. §. In 
both cases, the renormalization of the 3-d rota- 
tion ("rfi") term is performed perturbatively us- 
ing Q. For the D meson, only NP-tad is used 
since the approximations in || are not applicable 
for smaller masses. 

In the fat-link clover cases, csw is set equal to 
the tadpole improved tree level value 1/mq. The 
light-light renormalization coefficients are taken 
from the perturbative calculations of Ref. . The 
heavy-lights (for which perturbative calculations 
do not exist) are normalized using the static-light 
results of 0. This should be roughly correct for 
the large values of aM at the B meson. 

For Wilson valence quarks, we again use the 
EKM norm and identify as the physical mass. 
However since the magnetic and kinetic masses 
are not equal in this case there is little point 
in including the 3-dimensional rotations, and we 
do not. The associated systematic errors are es- 
timated as in Q. We normalize the heavy-light 
axial current as a function of mass perturbatively 
||. For the central values, the static-light value 
of the scale, (?g L « 1.43/a is used. (This scale 
is different from that quoted in Ref. g and used 
by us previously |Q; for further discussion see Q.) 

For the light quarks, 3 - 5 n values are used, 
in the range 0.7 - 2.0 m s . To get to physical 
light quark masses, we extrapolate vs. m-i 
quadratically and f^^m,Q q and /q 9 vs. m2 lin- 
early. The systematic error of the chiral extrap- 
olation is estimated by comparing with results 
from quadratic extrapolations for all the above. 
For heavy clover quarks, we use 3 - 5 k values 
(giving meson masses in the range 1.8-5 GeV); 
in the Wilson case we have ~ 10 heavy quark 
masses as well as static heavy quarks. Contin- 
uum extrapolation is done using both constant 
and linear fits, and the spread is taken as an es- 
timate of the systematic error. 

Figure [l] shows the continuum extrapolations 
of our full set of data for Jb, for both quenched 
and dynamical lattices. Both the NP-IOY and 
NP-tad quenched calculations, which should have 
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Figure 1. The complete set of Jb data, for both 
quenched and dynamical lattices, showing the dif- 
ferent continuum extrapolations. The line for fit 
(B), which is very close to the line for fit (A), has 
been shifted downward slightly for clarity. 



small discretization errors, are extrapolated with 
constants. The results are consistent with both 
the linear and constant quenched Wilson extrap- 
olations, although they favor the constant extrap- 
olation. All four extrapolations are averaged for 
our quenched central values; the spread deter- 
mines our discretization error. 

The dynamical lattice data clearly favor a con- 
stant fit; the best linear fit has a tiny slope. 
However, to be conservative, we also make a lin- 
ear extrapolation from the data for the smallest 
a, with the same slope as the linear fit for the 
quenched set. Other systematic errors, including 
the extrapolation in 1/Mq 9 , finite size errors, ef- 
fects of excited states, higher order perturbative 
corrections, etc., are estimated in basically the 
same way as in Refs. @^]. However the range 
of q* values considered for the perturbative error 
has changed along with the change in the central 
value of q*. We now consider 1/a < q* < 2q^ L , 
where q* SL is the value calculated in . 

As reported in Ref. the fat-link clover re- 
sults for decay constants Jd,d,.b.b s on dynamical 
lattices were found to be much smaller than the 
apparent continuum limit results of the Wilson 
quarks, thus raising the possibility that at least 



one of these results had a large uncontrolled sys- 
tematic error. To study this issue, we have com- 
puted with clover fermions on quenched lattices 
(set CP2) with four different levels of fattening: 
links fattened 2, 6 and 10 times with a fattening 
coefficient c=0.45, and links fattened 7 times with 
c = 0.25. A comparison with the thin- link clover 
computations (from sets CP1 and F) is shown in 
Fig. ||a. The fat-link results are considerably sup- 
pressed compared to those from the thin links. As 
mentioned above, the latter are consistent with 
the results for continuum-extrapolated quenched 
Wilson fermions. Furthermore, in an unpublished 
direct investigation of the static qq potential, we 
found that the form of the qq potential at short 
distances is distorted by fattening. Figure ||a 
also reveals the interesting fact that even the low- 
est level of fattening studied by us, viz., 2 levels 
of smearing with a fattening coefficient of 0.45, 
already substantially suppresses the decay con- 
stant. In fact, there is not much difference in the 
values of the decay constants between the four 
different levels of fattening we studied. 

The fat-link studies on quenched lattices, where 
the approach to the continuum limit is in much 
better control, allow us to attribute the cause of 
the discrepancy between Wilson and fat clover re- 
sults on dynamical lattices to fattening. We note 
however that the fattening effects in our calcu- 
lation include not only the change in the short- 
distance potential, but also the possible error 
in using the static-light rather than heavy-light 
renormalization. Be that as it may, the quenched 
studies clearly show that we should not average 
the dynamical fat clover and Wilson results (as we 
did in Ref. Q); instead, we use the Wilson val- 
ues only. Interestingly, correcting the fat-clover 
dynamical-lattice result by a factor obtained from 
the quenched lattices at comparable lattice spac- 
ing produces agreement with the Wilson fermion 
results (Fig. ^b), indicating that the suppression 
of decay constants due to fattening has similar 
origins for the dynamical and quenched lattices. 

Our final results for heavy-light decay con- 
stants in the quenched approximation are (in 
MeV for the decay constants): 



f B = 173±6± 16; 



/ B .= 199±5±*j; 
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Figure 2. a) fs for clover fermions on quenched 
lattices, b) fs on dynamical lattices, for both 
Wilson and fat-link clover fermions. Also shown 
is the "corrected" fat-link result. 

= 200±6±if; /d.= 223±5±#; 
J -p- = 1.16 ±0.01 ±0.02; 

IB 



f'D 



1.13 ±0.01 ±0.02. 



The errors are statistical and systematic (within 
the quenched approximation), respectively. The 
results differ from those in Ref. [Q due to: 

1) The inclusion of new data from sets CP, CP1 
and F. 

2) Setting the central value of the heavy-light 
scale as q* SL from jjj, rather than from ||. 

3) Other changes in analysis, motivated by the 
new runs. These include the error estimate for the 
chiral extrapolation (some alternative fits used 
previously are excluded by the new data) and 
the central value of the continuum extrapolation 
(averaging our four possible versions rather than 
taking only the linear Wilson fit). 

Our results for the dynamical lattices are (in 
MeV for decay constants): 



f D = 215 ± 5t£±g; f Ds = 241 ± 4±|?±g; 

f B 
fD. 
fD 



1.16 ±0.01 ±0.02 ±0.02; 



i.i4 ± o.oi±8:8| ± o.o2. 



Here the errors are statistical, systematic within 
Nf — 2, and systematic (due to partial quenching 
and missing the virtual strange quark), respec- 
tively. The last error is taken to be the largest 
of four components : 1) half the difference of the 
quenched and the Nf = 2 results (see Ref. [2)), 
2) the spread obtained from setting the scale by 
m p instead of by f n , 3) for quantities involv- 
ing the strange quark, the change in fixing k s 
from instead of the pseudoscalars, 4) the dif- 
ference from the "full unqucnching" result with 
TOvaionco = ^dynamical- For the individual decay 
constants, errors 1) and 2) are the largest; for the 
ratios, however, errors 3) and 4) dominate. 

REFERENCES 

1. C. Bernard et al, Phys. Rev. Lett. 81 (1998) 
4812. 

2. C. Bernard et al, Nucl. Phys. B (Proc. Suppl.) 
83-84 (2000) 289. 

3. A. El-Khadra et al, Phys. Rev. D 58 (1998) 
014506. 

4. M. Liischer, et al, Nucl. Phys. B491 (1997) 
323; ibid, 344. See also R.G. Edwards, U.M. 
Heller and T.R. Klasscn, Phys. Rev. Lett. 80 
(1998) 3448. The constant 6 a is the perturba- 
tive value (S. Sint and P. Wcisz, Nucl. Phys. 
B502 (1997) 251) with q* = 1/a (which re- 
produces the nonperturbative value for by)- 

5. K-I. Ishikawa, T. Onogi and B. Yamada, 
Nucl. Phys. B (Proc. Suppl.) 83-84 (2000) 
301, and private communication. 

6. C. Bernard, these proceedings. 

7. T. DeGrand and C. Bernard, Nucl. Phys. 
B (Proc. Suppl.) 83-84 (2000) 845, and in 
preparation. 

8. Y. Kuramashi, Phys. Rev. D 58 (1998) 
034507. 

9. O. Hernandez and B. Hill, Phys. Rev. D50 
(1994) 495. 



f B = 191 ± Gtlitl 1 ; f B . = 217 ± 5tgtS 



+33+9. 



